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ABSTRACT 

The discovery of quasars with redshifts higher than six has prompted a great deal of 
discussion in the literature regarding the role of quasars, both as sources of reioniza- 
tion, and as probes of the ionization state of the IGM. However the extreme ultra-violet 
(EUV) spectral index cannot be measured directly for high redshift quasars owing to 
absorption at frequencies above the Lyman limit, and as a result, studies of the impact 
of quasars on the intergalactic medium during reionization must assume a spectral en- 
ergy distribution in the extreme ultra-violet based on observations at low redshift, 
z < 1. In this paper we use regions of high Lya transmission (near-zones) around the 
highest redshift quasars to measure the quasar EUV spectral index at z ~ 6. We 
jointly fit the available observations for variation of near-zone size with both redshift 
and luminosity, and propose that the observed relation provides evidence for an EUV 
spectral index that varies with absolute magnitude in the high redshift quasar sam- 
ple, becoming softer at higher luminosity. Using a large suite of detailed numerical 
simulations we find that the typical value of spectral index for a luminous quasar at 
z 6 is constrained to be a = 1.3^03 for a specific luminosity of the form L„ oc v~ a . 
We find the scatter in spectral index among individual quasars to be in the range 
Aa ~ 0.75 — 1.25. These values are in agreement with direct observations at low red- 
shift, and indicate that there has been no significant evolution in the EUV spectral 
index of quasars over 90% of cosmic time. 

Key words: cosmology: theory - galaxies: high redshift - intergalactic medium - 
quasars 



1 INTRODUCTION 



The discovery of distant quasars has allowed detailed ab- 
sorption studies of the state of the high redshift intergalac- 
tic medium (IGM) at a time when the universe was less 
than a billion years old ( |Fan et al.|2006||Willott et al.|2007 
|2009[ ). Beyond z ~ 6 several of these quasars show a com- 
plete Gunn-Peterson trough in their spectra blueward of the 
Lya line ( White et al.|2003| . However, the spectra of these 
distant quasars also show enhanced Lya transmission in the 
region surrounding the quasar, implying the presence of ei- 
ther an HII region (Cen & Haiman 2000), or of a proximity 
zone. The interpretation of these transmission regions has 
been a matter of some debate. Different arguments in favour 
of a rapidly evolving IGM at z > 6 are based on the proper- 
ties of the putative HII regions inferred around the highest 
redshift quasars ||Wyithe fc Loeb||2004| I Wyithe et al.||2005 



Mesinger et al.|2004[|Kramer fc Haiman|2009|>. On the o ther 
hand, IBolton fc Haehnelt| ( |2007a| >, |Maselli et al.| ( |2007| > and 



high redshift quasar spectra blueward of the Lya line could 
also be produced by a classical proximity zone. In this case, 
the spectra provide no evidence for a rapidly evolving IGM. 
Importantly, these, and other analyses rely on an assumed 
value for the EUV spectral index to convert from the ob- 
served luminosity (redward of rest-frame Lya) to an ionizing 
flux, which is the quantity of importance for studies of reion- 
ization. As a further example, a value for the EUV spectral 
index must be assumed to estimate the quasar contribution 



to the reionization of hydrogen (e.g. Srbinovsky & Wyithe 
20071. 



The spectral energy distributions of luminous quasars 
are thought to show little evolution out to high redshift ( Fan 
|2006[ ). For example, broad emission line ratios at z ~ 5 have 



similar values to those observed at low redshift (Hamann & 



Ferland 1993 Dietrich et al. 20031. Moreover, optical and 



Lidz et al. (20071 have demonstrated that the features in 



infrared spectroscopy of some z ~ 6 quasars has indicated a 
lack of evolution in the optical-UV spectral properties red- 
ward of Lya ( |Pentericci et al. 2003 Vanden Berk et al.|2001[ ) . 
At higher energies, the optical/IR -to- X-ray flux ratios (e.g. 
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Brandt et al. 



2002; Strateva et al. |[2005 1 and X-ray contin- Here the near-zone size, i? CO r, has been corrected for lumi- 



uum shapes ( Vignali et al. 2003 1 show at most mild evo- 



lution from low redshift. However Jiang et al. (20101 have 



recently reported hot-dust-free quasars at z ~ 6 which have 
no counterparts in the more local Universe. These quasars 
are thought to be at an early evolutionary stage, and offer 
the first evidence that quasar properties have evolved since 
the reionization era. 

As yet, the spectral energy distribution blueward of the 
Lyman limit, which is critical for studies of reionization, is 
not measured at z ~ 6 since the quasar spectra are subject to 
complete absorption. In the future the EUV spectral index 
could be measured from 21cm observations of the thickness 
of ionizing fronts of the HII regions around high redshift 
quasars ( Kramer & Haiman|2008l) . However in the meantime 



the assumed value of the spectral index is usually based on 
direct observations of the EUV spectrum of quasars at z < 1, 
where they are not subject to absorption from the IGM. 

In this paper we show that the extent of Lya transmis- 
sion regions around high redshift quasars is sensitive to the 
EUV spectral index, and so provides a probe of this other- 
wise hidden portion of the spectrum. Our paper is presented 
in the following parts. We first summarise quasar near-zones 
(§ [2j , the observed near-zone relation (§[3|, and our associ- 
ated numerical modelling (§[4|. We then present constraints 
on the relation between near-zone size, quasar redshift, and 
quasar luminosity (§ p5j|. In § [6] and § [7] we present a dis- 
cussion of the astrophysical interpretation of our results. As 
part of our analysis we provide a revised estimate of the 
ionizing background at z = 6 and jointly constrain the ion- 
izing background and EUV spectral index. Some additional 
uncertainties in our modeling are discussed in § [8] and our 
conclusions are presented in § [9] In our numerical examples, 
we adopt the standard set of cosmological parameters (Ko- 
matsu et aL]|2009[ ), with values of Q h = 0.044, ft m = 0.24 
and Oa = 0.76 for the matter, baryon, and dark energy 
fractional density respectively, and h = 0.72, for the dimen- 
sionless Hubble constant. 



2 NEAR-ZONES 

In view of the difficulty of determining the physical condi- 
tions of the quasar environment, as well as the observational 
challenges associated with measuring the Lya transmission 



at high redshift, Fan et al. (20061 defined a specific radius, 
hereafter referred to as the near-zone radius, at which the 
normalised Lya transmission drops to 10 per cent after be- 
ing smoothed to a resolution of 20 A. This value of 10 per 
cent is arbitrary, and is chosen to be larger than the average 
Gunn- Peterson transmission in z ~ 6 quasar spectra (< 4 
per cent). In order to make the measurement unambigu- 
ous, the near-zone radius is defined at the point where the 
transmission first drops below this level, even if it rises back 
above 10% at larger radii. The measured near-zone radius is 
therefore dependent on the spectral resolution. 



Fan et al. (2006) found a striking relation between near- 



zone size and redshift, which they quantify using the expres- 
sion 



nosity according to the relation 

iW = R'27 X W -0M27+M li5o) / Bj 



(2) 



where B = 3 and R27 is the mean size of a near-zone at z — 6 
around a quasar of M1450 = —27. The value of B = 3 (which 
is roughly consistent with observations, Carilli et al.|2010 1 is 
motivated by the evolution of an HII region ( |Haiman|2002[ ) 
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where /hi is the fraction of hydrogen that is neutral, N is the 
rate of ionizing photons produced by the quasar, and t is the 
quasar ageQ The size of the near-zones is found to increase 
by a factor of 2 over the redshift range from 6.4 > z > 5.7, 
from which IFan et al. ( 2006 1 inferred the neutral fraction to 



increase by an order of magnitude over that time. 

On the other hand, if the near-zone size is set by res- 
onant absorption in an otherwise ionized IGM rather than 
by the boundary of an HII region, then its value is instead 
approximated by the expression 
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where Au m is the normalised baryon density in the trans- 
mitting regions as a fraction of the cosmic mean (rn m = 2.3, 
corresponding to the 10 per cent transmission), and T is the 
IGM temperature ( jBolton fc Haehnelt||2007b[). The param- 
eter a is the EUV spectral index. Equation Q differs from 
equation Q in three important ways. Firstly, if near-zone 
size is set by resonant absorption then it is not sensitive to 
the quasar lifetime. Secondly, resonant absorption results in 
a power-law index describing the dependence of near-zone 
size on luminosity with a value of B — 2 rather than B = 3. 
Thirdly, the near-zone size is sensitive to a, suggesting that 
near-zone size provides an opportunity to measure this un- 
known parameter. 

The precise definition of near-zone radius suggested by 



Fan et al. ( 2006 I allows for quantitative comparison with 



simulations. Wyithe et al. (20081 and Bolton et al. (20101 



-Rp = Rear + A x(z - 6). 



(1) 



have modelled near-zones using hydro-dynamical simula- 
tions with radiative transfer, combined with a semi-analytic 
model for the evolving, density dependent ionizing back- 
ground. The modelling in the latter work is able to produce 
both the amplitude and evolution of the near-zone sizes, as 
well as the statistics of absorption lines within the near- 
zone, and is based on quasars in a highly ionized IGM. 
In this regime it is equation Q with B — 2, rather than 
equation |3| with B = 3 that is applicable. This indicates 
that a varying ionizing background could be responsible for 
the observed relation, providing an alternative scenario to a 



1 This calculation assumes clumping and recombinations to be 
unimportant, that the density is at the cosmic mean over the large 
scales being considered, and that the quasar lifetime is much less 
than the age of the Universe. 
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Figure 1. Measured sizes _R P plotted against M1450 (Left) and 2 (Right) respectively, for observed (black dots) and simulated (triangles) 
near-zones. Measured near-zone sizes are also plotted for CFHQS J1509-1749, CFFQS J2329-0301 and CFHQS J0210-0456 (squares). 
The simulated data has been computed at the best fit in ojq and a\ (see text for details). 



strongly evolving neutral fraction. Unfortunately, this sce- 
nario implies that the evolution of quasar near-zones cannot 
necessarily be used to infer an evolution in neutral fraction. 
Thus, although intriguing, it could be argued that the ob- 
served near-zone relation has not yet taught us anything 
concrete about reionization. 

When discussing near-zones most authors have concen- 
trated on the relation between size and redshift, with mea- 
sured size corrected for luminosity dependence as part of 
this process. In most cases these studies used the scaling 
of near-zone in proportion to luminosity to the one third 
power (e.g. Fan et al.|2006 Carilli et al.|2010 |. As described 
above this value is appropriate if the near-zone is expanding 
into a neutral IGM. On the other hand, a power of one half 
(B = 2) is more appropriate if the near-zone corresponds 
to a radius where resonant absorption in an ionized IGM 
results in the 10 per cent transmission ( |Bolton fc Haehnelt| 
2007a I. Thus the correct scaling between near-zone radius 



and quasar luminosity is uncertain. In this work we there- 
fore treat the power-law index B as a free parameter in a 
2-dimensional relation describing near-zone radius as a func- 
tion of quasar luminosity and redshift. We find that current 
near-zone data places strong constraints on the value of B. 
Moreover, rather than focus on the interpretation of the pa- 
rameter A with respect to the end of hydrogen reionization, 
in this paper we concentrate instead on understanding the 
implications of the parameters B and R27 for the EUV spec- 
tral index blueward of the Lyman limit for z ~ 6 quasars. We 
show that the near-zone relation offers a probe of the EUV 
spectrum of quasars at these early epochs, information that 
is otherwise concealed by the Lyman-limit absorption, and 
so has previously only been studied for quasars at redshifts 
of z<\ |Telfer et al.]2002 |. 



3 OBSERVED NEAR-ZONE SAMPLE 



Carilli et al. ( |2010 1 have assembled a sample of 25 z ~ 6 
quasars which have quality rest-frame UV spectra and red- 



shift measurements (Fan et al. 2006 Jiang et al. 2008 1. Eight 
of the sample are detected ( |Wang et al.||2010 i in CO. An- 
other nine have redshifts measured ( Kurk et al.|2007 Jiang 
et al.|2007[ ) from the Mg II line emission. For the other eight 



objects, Carilli et al. (20101 adopt redshifts from the rele- 
vant discovery papers, which are mainly determined with 
the Lya+NV lines ( Fan et al.|2006 1. The measurements are 
summarized in Table 1 of Carilli et al. ( 2010| ), including the 
quasar absolute AB magnitude at 1450 A (M1450), quasar 
redshifts (z), and near-zone size (-R p ). Three sources in the 
original sample ( Wang et al.|2010] ) listed in Table 1 of |Carilli| 
eUuTlJlHOl, J0353 + 0104, J1044-0125, and J1048 + 4637, 



are broad absorption line quasars (Jiang et al. 2008 



et al.|2006 1, while the source J 1335 +3533 has a lineless (Fan 



Fan 



et al. 2006 1 UV spectrum. We exclude these quasars from 



our analysis owing to the fundamentally different nature of 
their intrinsic spectra ( Carilli et al.poTo Fan et aTf 2006). 
Thus the sample consists of 21 z ~ 6 quasars which we use 
to analyse the near-zone relation. Near-zone sizes are plot- 
ted in Figure[l]as a function of both absolute magnitude and 
redshift. Typical measurement errors ( Carilli et al.|2010 1 are 
1.2 Mpc, 0.4 Mpc, and 0.1 Mpc in the estimated near-zone 
size (-Rp) introduced by UV, Mgll, and CO-determined red- 
shift uncertainties, respectively. There are clear trends with 
both luminosity ( Carilli et al.| 2010) and redshift (Fan et al. 
2006} [Carilli et al.|2010"f ~ 

We note that quasar near-zone sizes have also been mea- 
sured ( |Willott et al.||2007| [Mof for CFHQS J1509 - 1749, 
CF FQS J2329 - 0301 and CFHQS J0210 - 0456. Follow- 
ing jCarjnTeFlir] (|20To| we do not include these near-zones, 
which are drawn from a different data set, in our analysis. 
However as may be seen from Figure [I] these near-zones lie 
on the Carilli et al.] ( 2010 ) correlation, and we have checked 



that their addition does not significantly alter our results. 



4 MODELLING OF NEAR-ZONES 

Our near-zone simulations combine a semi-analytical model 
for the evolving, density dependent photo-ionization rate 
in the biased regions surrounding quasars, with a radiative 
transfer implementation, and realistic density distributions 
drawn from a high resolution cosmological hydrodynamical 
simulation. These simulations are discussed in detail else- 
where ( Bolton et al.|2010 l, and their description is not re- 
produced here. However for completeness we show the evolu- 
tion of the background photo-ionization rate for our fiducial 
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Figure 2. The fiducial model for the ionizing background. Left: The evolution of the ionizing background | |Wyithe et"aj"1|2008| ) . The 
grey curve shows the background in the mean IGM, while the solid, dotted and dashed black curves show the mean and 1-sigma range 
of ionizing background in the biased regions within 5 Mpc of a 10 13 Mq quasar host halo. The data points show our estimate of the 
ionizing background at z = 5 and z = 6 (error bars represent the 68% range). The point at z = 4 has been reproduced from [Bolton fc| 
IHaehnclt (2007ajl. Right: The ionizing background as a function of radius in our fiducial model (for a quasar observed at z = 6). The 



ionizing background has been computed as a function of proper time along the trajectory of a photon emitted by the quasar, rather than 
at the proper time of the quasar. 



simulations in the left hand panel of Figure[2] The grey curve 
shows the photo-ionization rate in the mean IGM, while the 
solid, dotted and dashed black curves show the mean and 
1-sigma range of the photo-ionization rate in the biased re- 
gions within 5 Mpc of a 10 13 Mq halo. Importantly, we con- 
struct the absorption spectra using an ionizing background 
computed as a function of proper time along the trajectory 
of a photon emitted by the quasar, rather than at the proper 
time of the quasar. This effect is appropriate if considering 
spectra at the end of the reionization era, when the ionizing 
background can evolve significantly during the light travel 
time across a quasar near-zone. The right hand panel of Fig- 
ure [2] illustrates the background photo-ionization rate as a 
function of radius in our model (for a quasar observed at 
z = 6). 

For each of the model near-zone sets constructed, the 
sizes of the simulated near-zones were calculated for a sam- 
ple of synthetic quasar spectra with the same redshifts and 
quasar absolute magnitudes as the |Carilli et al.| ([2010) com- 
pilation. The synthetic spectra are designed to match the 
spectral resolution (R ~ 2500) and signal to noise of the 
observational data (S/N = 20 per 3.5 A pixel) as closely as 
possible. For the fiducial simulations an EUV spectral index 
of a = 1.5 (where dL/dv oc v~ a ) was assumed to calculate 
the corresponding ionizing photon rate (N) from the ob- 
served M1450. Additionally, a temperatu re at mean density 
of T = 236001 moo K has been measured ( |Bolton et al.|2010 1 
within the near-zone around SDSS J0818+1722. Our fiducial 
simulations are consistent with this constraint, and we fur- 
ther assume this as the best fit temperature when modelling 
near-zones throughout our analysis. This level of complexity 
is important as the near-zone size ( Fan et al.|2006| is mea- 
sured based on the first resolution element where the trans- 
mission in the spectrum drops below 10 per cent, meaning 
that the Lya forest must be simulated in detail in order for 
the modelling to be considered reliable. Examples of some 
of the synthetic quasar spectra used in our analysis are dis- 
played in Figure [3] 
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Figure 3. Example synthetic near-zone spectra. In both panels, 
the dotted curve shows the spectrum after being smoothed to a 
resolution of 20A in the observed frame while the vertical dashed 
line marks the measured near-zone size. This corresponds to the 
position where the smoothed flux first drops below a transmission 
of 10 per cent. Top: Quasar at 2 = 5.8918 with M 1450 = -26.09. 
Bottom: Quasar at 2 = 6.247 with M1450 = —26.67. 



5 NEAR-ZONE RELATION CONSTRAINTS 

As noted in the introduction, the appropriate value of B 
describing the power-law relation between near-zone size and 
luminosity is uncertain. Therefore, rather than correct for 
luminosity by assuming a fixed B, and fit equation |l]) to 
determine the evolution of near-zone size with redshift as 
has been done previously, in our paper we instead model the 
observed near-zone plane using free parameters to describe 
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Figure 4. Projections of x 2 f°r observed (black) and simulated (grey) near-zone sets assuming the fiducial evolving model. We have 
introduced intrinsic scatter (cr b s = 1.86 Mpc and <r s ; m = 1.27 Mpc, respectively) so that the reduced \ 2 is unity for the best fit model. 
The contours represent A% 2 = 1, 2, 3 and 5 relative to the corresponding minimum value of x 2 - 



both the evolution with redshift (A) and the evolution with 
absolute magnitude (B), 

R P = R27 x 10 -°-4(27+m 145o) / s +Ax(z~6). (5) 

In this section, we use equation (|5| to quantify constraints 
from the near-zone plane at z ~ 6. 

We find the parameter sets that describe the near-zone 
plane by fitting equation |5| to both the observed near- 
zone relation and to our simulated data sets. The black 
and grey contours in Figure [4] show the resulting projec- 
tions of x 2 onto the parameter spaces (A, B) and (^4, R27) 
for the observed and fiducial simulated near-zones, respec- 
tively. When calculating \ 2 we introduce intrinsic scatter 
into equation ([HJl so that the reduced \ 2 °f t ne best fit 
model is unity. We find a smaller intrinsic scatter (o- s im = 
1.27±0.04 Mpc) in the simulated data, indicating that sight- 
line variation is not responsible for all of the observed scatter 
(cr obs = 1.86 ±0.06). 

While the simulations reproduce the size evolution with 
redshift, the luminosity dependence is not well described, 
with B ~ 2 for the simulations (equation [4| but B ~ 3 for 
the observed sample. There are two explanations for this 
disagreement between the data and fiducial model which we 
now discuss in turn. 



6 NEAR-ZONE RELATION DETERMINED BY 
HII REGIONS 

Firstly, the observed value of B = 3 is expected if the 
quasars were surrounded by HII regions embedded in a neu- 



tral IGM ( Cen & Haiman 2000 1. Naively, this could be inter- 



preted as evidence for a significantly neutral IGM. To test 
this, we calculate the expected scatter in the near-zone rela- 
tion for comparison with observation. In contrast to the case 
for a proximity zone in an optically thin IGM, the size of 
an HII region depends on the age of the quasar at the time 
of observation. As a result, there is scatter in the observed 
size even for fixed redshift and quasar luminosity. We can 
estimate the scatter due to the age of the quasars by noting 
that if there is equal probability of observing a quasar at 
any age during its lifetime, then the probability distribution 
for the HII region radius at fixed luminosity L and redshift 



z, given a quasar lifetime t n 
dP 



3R P 

R3 



(6) 



where i? max is the size reached after t max . Given this dis- 
tribution, the mean is (R p ) = 3/4i? max , and the variance 
is erfj = (Rp) — {Rp} 2 = (^p) 2 /15. Observationally, since 
we have (Rp) ~ -R27 ~ 8 Mpc, the scatter in HII region 
radius owing to random quasar age is <jr ~ 2 Mpc. Thus 
the scatter in observed quasar age accounts for all of the 
observed scatter (a b s = 1.86 ± 0.06) in the near-zone sizes. 
However in the case of HII regions we would also expect 
additional scatter in the size owing to other quantities like 
the quasar lifetime and the spectral index (we return to this 
point below), as well as scatter due to inhomogeneities in 
the density field and ionization structure along the differ- 
ent quasar sight-lines. In the case of the latter, the sizes of 



HII regions near the end of reionization (Furlanetto et al 



20041 are thought to vary over the range 1-2 proper Mpc, 



which would lead to a component of scatter in addition to 
the ~ 2 Mpc expected from the quasar age. Thus, the total 
scatter in an HII region generated near-zone relation would 
be in excess of o\r ~ 3 Mpc. We therefore infer that although 
quasar HII regions would lead to a near-zone relation with 
the appropriate value of B ~ 3, the observed relation be- 
tween near-zone size and magnitude would be too tight in 
this case. This scatter based constraint would be alleviated 
if the near-zone sizes corresponded to resonant absorption 
within an HII region (the scatter would be smaller since the 
size is independent of lifetime in this case). However, in this 
alternative case where resonant absorption sets the near- 
zone size, our modelling suggests that B ~ 2 rather than 
B ~ 3 should be observed. As a result, HII regions cannot 
provide a viable physical explanation for the observational 
data (see also |Maselli et al. 12009] ). 



7 CONSTRAINING THE EUV SPECTRAL 
INDEX 

An alternative scenario to explain this trend is provided by 
an EUV spectral index which is luminosity dependent, so 
that at brighter absolute magnitudes the ionizing flux is 
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smaller. To quantify this point, we first note that the ob- 
served quantities are the near-zone size and the absolute UV 
magnitude. Given this UV magnitude, the ionization prop- 
erties of the surrounding medium are dependent on the ion- 
izing luminosity in photons per second (N) , and the spectral 
index blueward of the Lyman break. The relation between 
N and a is 



N 



hv 



v 



(7) 



where L„ L is the luminosity at the Lyman limit i^y, h is 
Planck's constant and i mlI is the ratio between the fre- 
quency where photons cease to contribute to ionization and 
the Lyman limit (we take £ max = 100, though this choice 
does not affect our results). To estimate the value of N for 



the quasars in the Carilli et al. (20101 sample, the observed 
M1450 has previously (Bolton et al. 20101 been combined 
with an assumed a spectral index of a = 1.5. Our approach 
to investigate whether a variable spectral index is able to ex- 
plain the discrepancy in B is therefore to modify the value of 
absolute magnitude that corresponds to the ionization rate 
at the near-zone radius measured in the simulations. This 
allows us to investigate the near-zone relations that result 
from a range of spectral indicies without re-calculating the 
radiative transfer simulations. Beginning with equation Q 
and equation Q and assuming a spectrum of the form 



L v (v) — Z/1450 



/ ^1050 \ 
V ^1450 / 



^1050 



V^Ly 



(8) 



we adjust the absolute magnitude corresponding to the 
model quasars (at fixed near-zone radius) for a variation 
in spectral index using 



M1450 
1 - :r 

X — 



v,»,..-^i(sr^ 



4.5 
max 



1 _ T ~(«+3) 
J- A max 



(1.71 - 0.86a + 0.33a 2 - 0.05a 3 ) " 



(9) 



where Mb, 1.5 is the observed absolute magnitude used in 
the fiducial model. Here we note that simulations performed 
with a = 0.5, a = 1.0, a = 1.5, a = 2.0, and a = 2.5 yield 
near-zone temperatures of T ~ 2.81 x 10 4 K, T ~ 2.34 x 10 4 K 
, T ~ 2.06 x 10 4 K , T ~ 1.89 x 10 4 K and T ~ 1.76 x 10 4 K 
respectively. These values are fit (to within 1%) by 



T = 35400 - 17700a + 6757. la 



lOOOa 3 



(10) 



This dependence contributes the last product in equa- 
tion jgj. 

To quantify the possible dependence of spectral index 
on luminosity we parametrize the value of a as a function 
of M1450 using 



a = Q - ai(M 1450 + 27), 



(11) 



where ao represents the mean spectral slope of quasars with 



M 



-27 and ot\ 



describes the variation of 



dMi450 

spectral index with magnitude. Note that since a is a func- 
tion of M1450, we evaluate it not at the observed M1450, but 
instead use an absolute magnitude M1450 that is adjusted 



to account for the average index ao at M1450 = —27, 

-1.8 

M' 1Am = M145C1.5 - 2.51og 10 



/ I>912 V 
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0.05ag) 



(12) 



We use these equations to constrain the likelihood of pa- 
rameter combinations (ao,ai) through convolution of the 
resulting joint probability distribution for the parameters 
A, B and R27 that describe simulated near-zone relations, 
with the corresponding relation for the observations 



where 
and 



dA dB dR 2 



d 3 P ob 



d 3 P„, 



d :i P ■ 

u 1 sin 



dAdBdR 27 dAdBdR 27 



Co 



dA dB dR 2 v 



d 3 P ■ 

u J sin 



(A,B,H 27 ) 



dAdBdR 27 " " ^ 13 ^ 

Here we have assumed flat prior probability distributions for 
each of the parameters A, B and R 2 j. As part of this proce- 
dure we calculate the intrinsic scatter for each set of param- 
eters (ao,ai), and impose the condition that this scatter be 
smaller than the observed scatter by weighting the calcu- 
lated likelihood with 



Co 



oc 



g 2(0.06^) 

1 otherwise. 



where cr obs < a s 



(14) 

Here the value of 0.06 dex represents the uncertainty in cr b s . 

Contours of likelihood for the values of ao and ai based 
on our fiducial model are plotted in the lower left panel of 
Figure [5] We find that the preferred value of ao is smaller 
than the assumed a = 1.5. In addition, since the fiducial 
simulations under-predict B (Figure ku, the value of ai is 
found to be greater than zero with high confidence, so that 
more luminous quasars have softer spectra. This is consistent 
with the results of |Strateva et al.| ( |2005[ ) , who found that 
brighter quasars are relatively less luminous in the X-ray, 
and is also expected from theoretical modeling of the quasar 
spectrum (Wa ndel fc Petrosian||1988[ ). 

The corresponding contours (grey curves) of A% 2 in the 
parameter spaces (A, B) and (A, P27) for the simulated near- 
zone relations computed using the most likely values of ao 
and ai are shown in the upper panels of Figure[5j along with 
contours for the observed near-zone sample (black curves). 
The best fit simulations reproduce the observed values of 
average size R 27 and index B. The parameter A, which 
is sensitive to the evolution of the ionizing background, is 
also consistent with observations ( |Fan et al.|2006| |Bolton fc| 
Haehnelt|20 07b). The near-zone sizes for the best fit model 
are plotted in Figure [l] (triangles). 

7.1 Intrinsic scatter in the spectral index 

As mentioned previously in section 5, the intrinsic scatter 
in the simulated near-zones (cr s im = 1-25 for the best fit 
model) is smaller than the scatter <r b s = 1-86 in the ob- 
served near-zone relation. Since the numerical simulations 
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12 





Figure 5. Upper panels: Projections of \ 2 for observed (black) and simulated (grey) near-zone sets for the best fit simulations (in cxq 
and «i) assuming the fiducial evolving ionizing background. We have introduced intrinsic scatter so that the reduced x 2 is unity for 
the best fit model. The contours represent A\ 2 = 1, 2, 3 and 5 relative to the corresponding minimum value of \ 2 . Lower left panel: 
Contours of likelihood for the values of ao and a\ (at 60% and 10% of the peak likelihood). Lower right panel: Contours of scatter in 
spectral index a. The point with error-bars shows the corresponding parameters measured | |Telfer et al. 2002) at low redshift. 



already include line-of-sight density variations in the IGM, 
the missing scatter 



AR P = 



-a 2 

sim 



(15) 



could be provided by scatter in the spectral index among 
individual quasars, which is known to be significant at lower 
redshift ( Telfer et al.|2002 l. Using the relation for the miss- 
ing scatter in R p we obtain the magnitude of scatter in the 
spectral index that is necessary to simulate the observed 
near-zone scatter 



Act 



da dM 

JmIlr 



AR = a 



B 

R27 



(16) 



Contours of this scatter are shown in the right hand panel 
of Figure [5] Values for scatter in spectral index ranging be- 
tween Aa ~ 0.75 and Aa ~ 1.25 correspond to the range of 
parameters ao and cti which lead to models for the near-zone 
relation that are in good agreement with the data. 



modelling of quasar near-zones assuming an ionizing back- 
ground that is independent of redshift. We first take the 
z = 6 value from our fiducial model; Figure [6] shows con- 
tours describing the resulting constraints (grey curves). The 
results are very similar to the fiducial case, with the excep- 
tion of A, which shows that the redshift evolution of near- 
zone size in the constant background model is smaller than 
observed. This is consistent with the previous inference that 
the trend of near-zone size with redshift is being driven by 



the rising intensity of the ionizing background ( Wyithe et al 
2008) at z ~ 6. We have also tested sensitivity to the ion- 
izing background amplitude by repeating our analysis for 
16 different values ranging between 1/100 and %/TO times 
the fiducial model. In the left panel of Figure [7] we show 
the resulting range for ao as a function of the background 
photoionization rate Tis = Thi/IO^ 12 s _1 . This figure illus- 
trates that results for a depend on the assumed value of 
Ti2. As a result, we next constrain T12 following the work 



of Bolton & Haehnelt ( 2007b I 



7.2 Constant ionizing background models 

Our simulations of quasar near-zones assume a semi-analytic 
model for the ionizing background (see Wyithe et al.|2008 |. 
It is therefore important to ask whether the results in this 
paper pertaining to the dependence of near-zone size on lu- 
minosity (i.e. the constraints on B), and the resulting con- 
clusions regarding the quasar spectral index a are sensitive 
to this model. To address this issue we have repeated our 



7.2.1 The Ionizing Background at z = 6 



Fan et al. ( 2006 ) present eight values of transmission T (with 
uncertainty 07-) measured in redshift intervals of Az = 0.1 
centred on redshifts in the range 5.9 < z < 6.25, and 14 
values centred in the redshift range 4.9 < z < 5.15. These 
values are listed in Table [T] From these values we calculate 
the mean value of transmission (T) = 0.004 ± 0.002 and 
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Figure 6. Upper panels: Projections of \ 2 for observed (black) and simulated (grey) near-zone sets for the best fit simulations (in ao 
and ai) assuming the non-evolving ionizing background. We have introduced intrinsic scatter so that the reduced x 2 is unity for the best 
fit model. The contours represent A\ 2 = 1, 2, 3 and 5 relative to the corresponding minimum value for x 2 - Lower left panel: Contours of 
likelihood for the values of ao an d ct\ (at 60% and 10% of the peak likelihood). Lower right panel: Contours of scatter in spectral index 
a. The point with error-bars shows the corresponding parameters measured ( |Telfer e t al. 2002) at low redshift. 
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Figure 7. Left Panel: The best-fit value for c*o as a function of Ti2 (solid curve), together with the 68% range (dashed curves). 
Central Panel: The joint probability distribution for Ti2 and «o- The prior probability for Ti2 is assumed to be flat in the logarithm 
(ciP pr i or /dri2 <x I/T12). Right Panel: The joint probability distribution for «o an d Q l- I n the central and right panels the contours 
represent levels that are at 60% and 10% of the peak likelihood. The final estimate for the EUV spectral index is «o = 1-3^q'3- 



(T) = 0.12 ± 0.01 respectively. To estimate this uncertainty 
we have calculated the standard error on the mean trans- 
mission via a bootstrap analysis which includes the uncer- 
tainties on individual parameters. The values of (T) = 0.004 
and (T) = 0.12 are consistent with the values quoted pre- 
viously ( jBolton &: H achnclt 2007a) at z — 6 and z = 5 re- 
spectively, although note that this study uses slightly wider 
redshift bins and does not use the independent data of 
|Songaila| ( |2004[ ). However our estimates of the uncertainty 



are smaller, particularly at 2 = 6 where the previously 
quoted error (Bolton & Hachnclt 2007a) is a factor of ~ 2 
larger. This difference can be attributed to the fact that 
rather than estimate the standard error on the mean trans- 



Bolton & Haehnelt (2007a) used an inter-quartile 



range as an estimate for the error on the transmission values 
along various lines of sight. 

To estimate the corresponding background photo- 
ionization rate we next combine previously pub- 
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4.92 


0.1276 


0.0011 


4.95 


0.1139 


0.0022 


4.98 


0.1002 


0.0050 


5.01 


0.1268 


0.0117 


5.04 


0.1567 


0.0073 


O.UD 


U.l / DO 


n m on 


5.06 


0.1285 


0.0042 


5.06 


0.1509 


0.0224 


5.07 


0.0751 


0.0011 


5.08 


0.0530 


0.0025 


5.10 


0.0898 


0.0020 


5.11 


0.1650 


0.0030 


5.13 


0.1293 


0.0064 


5.13 


0.1243 


0.0050 


5.90 


0.0108 


0.0033 


5.93 


0.0125 


0.0022 


5.95 


0.0038 


0.0005 


5.95 


0.0060 


0.0010 


6.08 


-0.0071 


0.0020 


6.10 


0.0012 


0.0010 


6.10 


0.0051 


0.0005 


6.25 


0.0015 


0.0005 



Table 1. Table of transmissions T and uncertainty a-f in Az 
0.1 regions from |Fan et al.| j2006| . 



lished ( Bolton & Haehnclt 2007a I data with a Monte-Carlo 
method. We first determine the resulting distribution of 
effective optical depths from our transmission estimates, 



finding values of r e g = 5.5lo'.4 an d r e a = 2.1 



at z 



6 



and 2 = 5 respectively. Here and below we have quoted 
the median and the 68% range. We note that the value of 
T e ff is not defined in the small number of cases where the 
Monte-Carlo trial value of T is negative. In such cases we 
set T = 10 -6 and note that our results are insensitive to 
the exact value provided it is well below the observational 
limits in the Gunn-Peterson troughs. Given the distribution 
of effective optical depths, we determine the corresponding 
probability distribution for the normalized photoionization 
rate Fi2 = Thi/IO -12 s _1 . To estimate the uncertainties we 
have utilized scaling relations derived from hydrodynamical 
simulations (see Table 4 of Bolton fc Haehnelt| |2007a[ ) 
for Ti2, as functions of r e ff, the IGM temperature [which 
we take to be log 10 (T/10 4 K) = ± 0.3], and the slope of 
the temperature density relation. We also utilize scaling 
relations as a function of the cosmological parameters. 

At z = 5 we find T 12 = 0A7tn l, while at z = 6 we ob- 



tain Ti2 



0.18 



+0.18 



We note that the smaller uncertainty 



in r c ff relative to the previous estimate ( Bolton & Haehnelt 
2007a| allows us to quote a probability density for the mea- 
sured ionizing background (with a non-zero value preferred 
at 2a) rather than an upper limit. The revised estimates of 
the ionizing background at z = 5 and z = 6 are presented 
in Figure [2] 



bility distribution for Fi2 and qq (central panel of Figure [7]). 
Similarly, we obtain the joint probability for ceo and «i (right 
hand panel of Figure |7|. This procedure yields an estimate 
of oio = I.3I0 3. This constraint represents the primary mea- 
surement of our paper. 



7.3 Comparison with spectral index at low 
redshift 



Telfer et al. 



( 2002 1 measured the UV spectral index for 
2. By analysing the composite spec- 
1.76 ± 0.12 



Telfer 



quasars at z ~ 1 

trum they found a mean value of (a) 

al. ( 2002 1 also measured a for a sub-sample of indi- 
vidual quasars, allowing them to estimate the variation 
Qi = 0.3 ±0.1 with luminosity as well as the intrinsic scat- 
ter Aa = 0.76. These estimates are plotted alongside our 
model constraints in Figure [B"|8| We find that the parameter 
combination ao = 1.76 ±0.1 2 and Qi=0.3±0.1is consistent 
with our inference for these parameters at z ~ 6. Moreover, 
we find that our inferred scatter in a at z ~ 6 agrees well 
with the Aa = 0.76 found bylTelfer et all (120021). Thus, as- 



suming our fiducial model of the ionizing background, there 
is no evidence for evolution in the EUV slope of quasars with 
redshift out to z ~ 6. 



8 UNCERTAINTIES 

In addition to those already discussed, we have checked sev- 
eral other sources of uncertainty. Specifically, we have recal- 
culated our constraints including the assumption of a fidu- 
cial model with a — 0.5 rather than a — 1.5, a model where 
the quasar continuum has been placed systematically high 
by 5% on the data, and a model that includes an approx- 
imate treatment of Lyman limit systems. Each of these ef- 
fects is found to have a relatively small effect on the derived 
parameters for the spectral index relative to the quoted un- 
certainty, as detailed below. 



8.1 Near-zone temperature 

There is an uncertainty in the absolute temperature of the 
IGM within the quasar near-zone, as well as scatter in the 
temperature from one near-zone to another. These uncer- 
tainties will impact on our constraints for the EUV spectral 
index. Firstly, any scatter in temperature among different 
near-zones will contribute to scatter in the observed near- 
zone sizes. However this scatter will be degenerate with scat- 
ter in the spectral index (Aa). Thus the potential presence 
of scatter in temperature does not introduce uncertainty into 
our estimated mean value of Qo or cti- However, it does imply 
that our quoted Aa should be considered an upper limit. 

On the other hand, uncertainty in the absolute value 
of the temperature within the near-zones could introduce 
a more serious error. The temperature of the IGM in the 
quasar near-zone has two independent contributions which 



7.2.2 Joint Probability for V 12 and ao 

We use the constraints for Ti2, together with the likelihoods 
for ao and a± as a function of Fi2, to generate a joint proba- 



2 Note that the evolution of a with luminosity is dominated by 
radio loud quasars, which formed half the total sample of quasars 
for which individual measurements of a were made by |Telfer et al.| 
(2002 k 
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Figure 8. Figure showing the systematic uncertainties on the joint probability distribution for c*o an< f oi\. Left panel: The dark contours 
show constraints that have been computed based on a model with a = 0.5. Central panel: The dark contours show the constraints in a 
case where the continuum has been deliberately placed at a level that is 5% too low. Right panel: The dark contours show the constraints 
in a case where where self-shielded absorbers have been added to regions with density contrast A > 16. In each case the model assumes 
an unevolving ionizing background, for which the light contours in each panel show the constraints. The contours represent levels that 
are at 60% and 10% of the peak likelihood respectively. 



add linearly, and hence the near-zone temperature has two 
components of uncertainty. Firstly, the mean IGM has an 
uncertainty in temperature resulting from heating during 
hydrogen reionization. This component of uncertainty in 
the temperature due to hydrogen reionization is already in- 
cluded ( jBolton fc Ha chnclt 2007a) in the calculation of the 
ionizing background. 



However in addition the reionization of Hell by the 
quasar may heat the IGM in the near-zone to well above 



this mean IGM temperature (Bolton et al. 20091. In our 



modelling we have assumed a temperature (Bolton et al. 
|2010[ ) that is consistent with observations of a single z ~ 6 
quasar (SDSS J0818+1722), and that the quasar is respon- 
sible for reionization of Hell within the near-zone. This lat- 
ter assumption means the temperature increase from Hell 
reionization by the quasar is directly related to the spectral 
hardness of the quasar, and hence to a, which is the quan- 
tity we are trying to measure. As a result, the assumption 
of an incorrect temperature could bias the inferred value of 
a in our analysis. A hardening of the spectral index leads 
both to higher temperatures within the near-zone and to an 
enhanced photo-ionization rate. By assuming a fixed tem- 
perature for all modelling our analysis would include only 
the latter enhancement, whereas both ionization and tem- 
perature influence the optical depth. We have therefore ac- 
counted for this effect via the scaling in equation (JoJ based 
on a calibration of near-zone temperature with spectral in- 
dex estimated from additional radiative transfer simulations. 
Furthermore, to check the importance of the temperature ef- 
fects, and at the same time check that our assumed fiducial 
value of a does not influence the results, we have repeated 
the constraints on ao and ot\ based on the fiducial unevolv- 
ing ionizing background using a simulated quasar with a 
value of a = 0.5 (and correspondingly hotter temperature). 
As shown in the left panel of Figure [8] we find that the as- 
sumed fiducial a (and corresponding fiducial temperature) 
does not strongly influence the constraints. 



8.2 Continuum placement 

In addition to errors in the redshifts of the quasars, uncer- 
tainty in continuum placement may also lead to an error in 
the near-zone radius R p . This error corresponds to an uncer- 
tainty ( Fan et al. 2006 I of At ~ 0.05, or equivalently to a rel- 
ative uncertainty in the transmission of AT/T ~ 5%. At the 
near-zone radius T = 10%, so that the uncertainty in trans- 
mission is AT ~ 0.005. Inspection of Figure 3 in |Carilli et al.| 
(2010) shows that when averaged over the near-zone sample, 
the gradient of transmission with radius is dT/dRp ~ 0.05 
Mpc - . Hence the uncertainty in the near-zone radius that 
is introduced by the uncertainty in continuum placement is 
AR P ~ dRp/dTAT ~ OTMpc. Thus, continuum placement 
is sub-dominant in terms of contribution to the uncertainty 
in near-zone size. We note that the uncertainty could be 
larger for individual quasars. However this would introduce 
additional scatter into the near-zone relation which would 
be degenerate with the other sources of scatter. To check 
this conclusion we have again repeated the constraints on 
«o and cvi based on an unevolving ionizing background us- 
ing simulated quasars for which the continuum has been 
systematically placed 5 per cent too low. As shown in the 
central panel of Figure [8] we find that the continuum place- 
ment has a systematic effect on near-zone size and hence 
Qo. However the shift is much smaller than the statistical 
uncertainty. We thus conclude that our results are robust to 
a systematic uncertainty in the continuum at this level. 



8.3 Lyman limit systems 

Finally, a potential shortcoming of our modelling is the as- 
sumption that the IGM is optically thin, which results in the 
absence of absorbers that are self-shielded with respect to 
the ionizing background. To address this we have tested the 
sensitivity of near-zone size by adding a population of self- 
shielded Ly-limit systems. We achieve this by setting all gas 
with density contrast A > 16 (corresponding to the density 
contrast of a clump that will self shield if immersed in our 



fiducial ionizing background, Furlanetto & Oh 2005) to be 
neutral prior to the quasar turning on (Bolton & Haehnelt 



2007a). We find that their inclusion does not alter the pre- 
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dieted distribution of near-zone sizes in most cases, and that 
as shown in Figure |J the presence of self-shielded systems 
should not significantly modify our measurement of ao rel- 
ative to other uncertainties. 

A related question concerns an alternative explanation 
for the observed B = 3 due to an enhanced number density 
of absorbing systems in the biased environments surround- 
ing luminous quasars. Such an enhancement could preferen- 
tially retard the observed near-zone sizes, leading to an in- 
crease in near-zone size with luminosity that is slower than 
expected in the optically thin regime. However we have pre- 
viously (Bolton & Haehnelt 2007a) investigated the depen- 
dence of near-zone size on quasar host-halo mass, finding a 
negligible effect. This is because the absorbers that set the 
near-zone size are located between 30 and 100 co-moving 
Mpc from the quasar. Thus, the environmental dependence 
of absorber abundance is unlikely to explain the observed 
B = 3. 



9 CONCLUSION 

The discovery of a strong correlation of the sizes of near- 
zones around high redshift quasars and their redshift in the 
range 5.7<z<6A ( Fan et al.||2006| [Carilli et al.||2010| has 
prompted a series or studies aimed at understanding its im- 



plications for the tail end of reionization (Fan et al 



2006 



Bolton fc Haehnelt|200"7al |Wyithe et al.||2008| |Bolton et at] 



2010| |Carilli et al.||2010p ln this paper we have noted that 
the amplitude of the near-zone radius relation is sensitive to 
the assumed EUV spectral index of the quasar. We there- 
fore performed a large suite of numerical simulations of the 
near-zone spectra using radiative transfer through a hydro- 
dynamical model of the IGM in the quasar environment, 
and find that the observed near-zone relation can be used 
to constrain the UV spectral index of z ~ 6 quasars. We 
find that the typical value of spectral index for luminous 
quasars at z ~ 6 is a — l.Stg'^ (for a specific luminosity of 
the form L„ oc v~ a ), where the uncertainty includes both 
the statistical uncertainty in the near-zone relation as well 
as uncertainty in the ionizing background. Based on com- 
parison of the scatter in our simulations and the scatter in 
the observed near-zone relation we infer that the scatter in 
the spectral index is Aa ~ 0.75 — 1.25. 



Carilli et al. (20101 noted that the data are consistent 



with near-zone sizes that scale with quasar luminosity as 



Rr, 



L 1 / 3 , where B = 3, and that the scatter in the 



near-zone radius-redshift relation is reduced if this scaling 
is applied. However as pointed out by Bolton fc Haehnelt | 
(2007a I, the value of the power-law index B describing the 



scaling of near- zone radius with luminosity that is physically 
appropriate depends on whether near-zone sizes are set by 
HII region boundaries (B = 3) or resonant absorption in 
an ionized IGM (B = 2). Motivated by the possibility of us- 
ing this second parameter in the near-zone relation to gain a 
better understanding of its meaning, in this paper we jointly 
fit the available observations for variation of near-zone size 
with both redshift and luminosity. From this analysis we 
find that the observations have an index of B ~ 3, with 
an uncertainty that excludes a value of B = 2 with high 
confidence. In contrast to the observations, we find that our 
fiducial numerical simulations predict a value of B — 2. 



We have discussed two possibilities for this difference. 
Firstly our simulations are conducted in a highly ionized 
IGM. It is possible that although these simulations are able 



to describe the size and evolution of the near-zones ( Bolton 
et al.|2010 l, they do not represent physical conditions of the 
IGM at z ~ 6, which could instead contain hydrogen that 
is significantly neutral as previously argued in |Wyithe ~Ez\ 
Loebl (120041 andlWyithe et al.l d2005b. Although this would 



lead to the observed B = 3, we do not believe that this is 
the case because the observed relation is too tight to be ex- 
plained by HII regions, the sizes of which contain a range 
of additional sources of scatter. Instead, we argue that the 
observed value of B ~ 3 is evidence for a UV spectral in- 
dex that varies with absolute magnitude in the high redshift 
quasar sample, becoming softer at higher luminosity. This 
finding is in agreement with theoretical modeling of quasar 
spectra ( Wandel fc Petrosian|1988 1. 

The results of this paper provide the first constraints 
on the properties of the EUV spectral index of the highest 
redshift quasars, which are not observable owing to absorp- 
tion blueward of the Lyman limit. The values we find are 
in agreement with direct observations at low redshift, and 
therefore indicate that there is no evidence for evolution in 
the EUV properties of quasars over most of cosmic time. The 
new constraints will aid all studies of reionization involving 
high redshift quasars, which have previously relied on an as- 
sumed value of the mean EUV spectral index measured at 
lower redshift. 
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